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Abstract. We consider lattice gas diffusive dynamics with creation- annihilation in the 
bulk and maintained out of equilibrium by two reservoirs at the boundaries. This stochastic 
particle system can be viewed as a toy model for granular gases where the energy is injected 
at the boundary and dissipated in the bulk. The large deviation functional for the parti- 
cle currents flowing through the system is computed and some physical consequences are 
discussed: the mechanism for local current fluctuations, dynamical phase transitions, the 
fluctuation-relation. 

1. Introduction 

In different contexts like biological systems or chemical reactions, many out of equilibrium 
dynamics combine diffusive and dissipative features which can be described macroscopically 
by reaction-diffusion equations. In order to analyze these systems from a microscopic point 
of view, a lattice gas model with fast stirring and a local creation-annihilation mechanism 
has been introduced in |DFL] by combining a symmetric simple exclusion process (SSEP) 
and a Glauber dynamics. Another example of driven dissipative system is a granular gas 
where the energy is dissipated through inelastic collisions. Direct analytic approaches of 
the granular gases are notoriously difficult, but stochastic toy models have been proposed in 
order to capture some aspects of the physics [HI ILLl ISLl [Fl IFPj or for simulation purposes 
|GRWj . The particle system introduced in |DFL] in contact with two heat bath shares also 
similar features. 

In non-equilibrium statistical mechanics, there is no analogous to the Gibbs formalism and 
the large deviation functional can be viewed as a substitute for the free energy. In particular, 
the large deviations of the current received a lot of attention over the last few years (see 
|BDGJL2l [D| [G] for reviews). Exact expressions of the large deviation functional have been 
obtained for conservative diffusive stochastic dynamics in contact with reservoirs [BDGJLl] 
|BDLj . In this paper, we generalize these results to the non-conservative stochastic dynamics 
which combine a fast stirring and a Glauber dynamics as in [DFLj . The relevant macroscopic 
parameters to describe these stochastic systems are the density, the conservative current 
(from the stirring dynamics) and the non-conservative current (from the Glauber dynamics). 
The large deviation functional (13.11) gives the exponential cost of observing a deviation of 
these three parameters simultanously. The large deviations for the density have been already 
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studied in |JLVj and tlie long-range correlations of the steady state in |BJj . The large 
deviations of the injected power were computed for a different non-conservative stochastic 
dynamics in |FPj . 

The large deviation functional derived in (13.11) measures the cost of observing a joint 
deviation of the density and the currents. Thus the density large deviations obtained in |JLVj 
can be recovered from the functional (13.21) by optimizing the deviations on the currents. In 
the same way, observing an atypical conservative current imposes a non trivial constraint 
on the density and the non-conservative current. As a consequence, we argue in section HTT] 
that the local conservative current fluctuations are produced by a different mechanism than 
in the case of conservative dynamics. 

For some choice of the Glauber rates, the steady state may exhibit a phase transition. 
This instability of the steady state leads to a new kind of dynamical phase transition for the 
current which is presented in section 14. 2[ 

The fluctuation relation for the entropy production has been investigated in the framework 
of dissipative systems in |BGGZ[ IFF] . Using the exact expression of the large deviation 
functional (13. 2p . we discuss its symmetry properties with respect to time reversal and recover 
the fluctuation relation for the entropy production. Contrary to the conservative dynamics, 
the entropy production is no longer proportional to the current injected in the system, but 
depends in a complex way on the currents and the density. 

2. Models and Notation 

To fix the idea, we first recall in section [TT] the microscopic model introduced in [DFLj . As 
it will be clear from section fI72i the details of the microscopic evolution will not be important 
to study the macroscopic limit and one can consider more general dynamics than the one 
introduced in section 12.11 The relevant macroscopic parameters are defined in section 12.21 
and some concrete examples are then discussed in section 12.31 

2.1. A microscopic model. We consider a particle system on a one dimensional chain 
with sites {1, A"}. The state of each site i can be occupied or empty and is encoded by 
rji G {0, 1}. During the infinitesimal time interval dt, each particle has a probability N'^dt of 
jumping to the left if the left neighboring site is empty, N'^dt of jumping to the right if the 
right neighboring site is empty. Furthermore at each site i, a creation or an annihilation can 
occur with probability c{i,ri)dt, where the rate c{i,r]) can depend on the local configuration 
around the site i. The simplest rates are constant creation and annihilation rates c~^,c~ in 
which case c{i,ri) = c'^{l — rji) + c'rji. We will also consider more complicated rates leading 
to dynamical phase transitions (I2.18p . Finally, to model the effect of reservoirs of particles 
at the boundaries another creation/annihilation mechanism is acting on the boundary sites: 
during each time interval dt a particle is created at site 1 (if it is empty) with probability 
N'^afdt and removed with probability N'^a^dt (if it is occupied). Similarly at site A^, 
particles are created or removed at rate N'^a^dt. 

We stress the fact that the stirring of the particles and the action of the reservoirs is 
much faster (by a factor A^^) than the creation/annihilation mechanism in the bulk. The 
particles perform random walks thanks to the stirring mechanism. Thus in the large A^ 
limit, a particle created at the boundary has a positive probability to go through the system 
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without being annihilated. This is the correct scahng to ensure a competition between the 
two mechanisms of the dynamics and therefore a non trivial thermodynamic limit (A^ — oo). 

Rephrased in mathematical terms, the previous dynamics is a Markov chain with generator 
given by 



where the configuration r/*'*+^ is obtained from 77 by swapping the occupation numbers of the 
sites + l and the configuration 77* is obtained from t] by changing the occupation number 
at site i from rji to 1 — rji. 

The microscopic currents associated to the previous dynamics are defined as: 
The conservative current. For a given edge + 1), let QT,{j,i+i) be the number of jumps 
from site i to site i + 1 minus the number of jumps from site z + 1 to site i during the time 
interval [0, T]. 

The non-conservative current. We denote by K^^ the number of particles created at site 
i during the time interval [0,T]. In the same way, K^- denotes the number of particles 
annihilated at i. 

Microscopically this leads to the following conservation law at each site i and any time t 



2.2. The hydrodynamic limit. One can also consider more general stirring processes than 
the SSEP of section [2TT] (e.g. the Zero range process or a Kawasaki dynamics [S]). In the 
macroscopic limit {N — ^ 00) this will lead to a non-trivial hydrodynamic limit provided that 
the diffusive dynamics is speeded up by a factor A^^ compared to the Glauber part. The 
details of the microscopic evolution will be averaged out in the thermodynamic limit and the 
macroscopic equations will depend only on a few relevant macroscopic parameters which we 
describe below. 

We first consider the purely diffusive evolution in contact with the reservoirs at the bound- 
ary (i.e. when c{i,r]) = 0). According to the hydrodynamic formalism [S], a given diffusive 
lattice gas can be characterized by the diffusion coefficient D{p) and the conductivity a{p) 
both depending on the density p. Following |BD1] , one way to define these coefficients is to 
consider a one dimensional system of length A^ connected to reservoirs at its two ends. Since 
the number of particles is preserved by the dynamics, it is equivalent to measure the current 
Qr,(i,i+i) through any bond {i,i + 1). When both reservoirs are at the same density p, the 
variance of the current Qr,(j,i+i) transfered during a long time T from one reservoir to the 
other is given for large A^ by 




Vt,i - Vo,i = i - K^ - + Qt,(i-i,i) - Qt,{i,i+i) ■ 



(2.1) 



lim 



{{QT,{i,i+l)) ) 

T 



a{p)N, 



(2.2) 
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where (■) denotes the expectation of the dynamics starting from the invariant measure. The 
conductivity (j{p) is defined by (12.21) . Note that the scahng wrt comes from the fact that 
the diffusive dynamics has been speeded up with a rate N'^. On the other hand if the left 
reservoir is at density p + 5 and the right reservoir at density p, the average current is given, 
for small 5 and large A^, by 

Um i^^|^ = D(p)iV5, (2.3) 

which is simply Pick's law and defines the function D{p). For the SSEP, a{p) = 2p(l — p) 
and D{p) = 1 [S]. 

At time t and for x in [0, 1], the macroscopic density is denoted by p{x, t) when N oo 

\i-Nx\ ^ £ 

We stress the fact that the microscopic diffusion has been speeded up by a factor A^^ thus 
t stands already for a macroscopic time. At the boundary, the density is imposed by the 
reservoirs at any time t 

p{M=Pa= p(l,t)=p,= _ , (2.5) 

where the reservoir densities are those of the dynamics defined in section 12. 1[ 
Finally we define the macroscopic conservative current at x G [0, 1] up to time t 

S <5i,(M+i)' (2-6) 

\i-Nx\ ^ e 

where e vanishes in the thermodynamic limit (A^ oo). Its time derivative Q{x,t) stands 
for the local current at x and according to (12. 3p . its typical value is given by 

Qix, t) = -D{p{x, t))d,pix, t) . (2.7) 
Thus the hydrodynamic limit is given by [S] 

dtp{x, t) = {D{p{x, t))d,p{x, t)) , (2.8) 
with the boundary conditions (12.51) . 

We investigate now the hydrodynamic behavior of the full dynamics. Let pp be the in- 
variant measure of the stirring dynamics at the uniform density p. (For the SSEP then pp 
is simply the Bernoulli product measure at density p). At time t = 0, the system starts at 
constant density p from the measure pp, then for any site i = [Nx] in the bulk 

dtpix,t = 0) = dMt = 0)) = -A{p) + C{p) , (2.9) 

where the macroscopic annihilation rate A and the macroscopic creation rate C are obtained 
by averaging the microscopic rates over the measure pp 

A{p)= pp{cii,r])7]i), Cip) = Pp{cii,7]){l-T]i)). (2.10) 
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The measure iip is in general no longer invariant for the full dynamics, however the fast 
stirring (with an A^^ rescaling) dominates locally and maintains the system in a local equi- 
librium. Even though, the density p(x, t) evolves under the action of the dynamics, the local 
statistics remain given by ^Jip{x,t) up to some small corrections which are vanishing for large 
A^. Thus when N oo, (12. 9p can be generalized to later times by adding the diffusion term 
(12. Sp and the density relaxes by following the hydrodynamic equation 

dtp{x,t) = d.,{D{p{x,t))dxp{x,t)) - V'{p{x,t)) , (2.11) 

where the potential is given by 

V'{u) = A{u) -C{u) , (2.12) 

and the boundary conditions are imposed by the left and right reservoirs p{0,t) = pa and 
p{l,t) = pb (12.51) . In particular, the steady state density p satisfies 

d,{D{p{x))d,p{x)) - V'{p{x)) =0. (2.13) 

These heuristic ideas have been rigorously justified in [DFLl IJLV^ IBLj for the microscopic 
model of section [2?T] (in this case D = 1). 

The macroscopic non-conservative current at x G [0, 1] up to time t is defined by 

\i-Nx\ ^ £ 

Its time derivative K{x, t) stands for the local current at x and is typically given by 

k{x, t) = -A{p{x, t)) + C(p(x, t)) . (2.15) 

After rescaling, the microscopic conservation law (12.11) reads 

dtp{x, t) = -dM^, t) + k{x, t) . (2. 16) 

2.3. Some concrete examples. The previous models can be used to describe a wide variety 
of physical systems. We describe below a few applications. 

Granular gasses. 

When there is no creation in the bulk, the dynamics introduced in section 12.21 can be 
viewed as a toy model to mimick the granular gases where energy is injected at the boundary 
and dissipated in the bulk by inelastic collisions. Similar dynamics have been considered 
in [HI ILLl ISLj . More precisely, for the model of section 12.11 if c{i,ri) = arji, then the 
hydrodynamic limit equation (12.111) reads 

dtp{x, t) = dlp{x, t) — ap{x, t) , (2-17) 

with boundary conditions Pa, Pb (12.51) . p{x,t) can be interpreted as a local energy dissipated 
at rate ap in the bulk and the stationary state (12.131) satisfies 

-f \ ^ \ , pb- paCosh{y/a) 

pix) = Pa cosh(Vax) H , . smh[y/ax] . 

smh(va) 



React ion-difFusion . 
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The microscopic model of section 12.11 was originally introduced |DFLj to model diffusion 
reaction type equations often used to describe chemically reacting systems. For a Glauber 
dynamics which satisfies detailed balance with respect to the Gibbs measure of a one dimen- 
sional Ising model with nearest neighbor interaction at inverse temperature /5, then 

c{i, r/) = 1 + 7(1 - 2r/,,)2(r/,_i + r],+^ - 1) + (-1 + 2T^i_^){-l + 2t^,+^)^^ , (2.18) 

with 7 = tanh(/5). The hydrodynamic limit (12. lip reads 

dtp{x,t) = dlp{x,t)-V'{p{x,t)), 

with the potential fl212|l 

V{p) = -(1 - 7)'(1 - P)P + 2f(l - P)V . (2.19) 

For 7 < 1/2, the potential V has a unique minimum, instead ^ is a double well when 7 > 1/2. 
In the latter case a phase transition may occur as there might be several solutions of fl2.13p 
which can be stable or not with respect to the fluctuating hydrodynamic. Depending on the 
value of the parameters and of the boundary conditions, the steady can be concentrated on 
more than one profile, leading therefore to a phase transition. Some consequences of these 
transitions on the current large deviations will be examined in section 14.21 

3. Large deviations 

3.1. The functional. We turn now to the large deviations of the currents and density. The 
corresponding functional can be completely characterized at the macroscopic level in terms 
of the coefficients a, D fl2.2ll2.3|l and A, C fIXTU]) . 

On average, the density follows the hydrodynamic equation (12. lip and the macroscopic 
currents {Q{x,t), K{x,t)} obey (12. 7p . (I2.15p . The probability of observing the evolution of 
an atypical density profile p{x,t) and atypical macroscopic currents {Q{x,t), K{x,t)} for 
< t < T is given for large by 



P[o,T] ({p,g,K}) ~exp [-NI[o,T]{p,Q,K)] (3.1 



where I[o,t] = cxd if the evolution does not satisfies the conservation law (I2.16P and the 
boundary conditions (12. 5p . Otherwise, X[o,t] is defined by 

,2 



. • \ (Q(x,t) + D(p(x,t))d^p(x,t)Y / N 

X[o,T](p,Q,if) = dt dx I ' ^ 2a{p{xt)) +^[p{x,t),Kix,t)^ 



(3.2) 



with 



'k^ + 4A(p)c(p) + k 



^p,K) = C{p) + A{p)-^JK^ + AA{p)C{p) + K\og\^^ | . (3.3) 

If C(p) = 0, then $ becomes 

00, if i: > . 
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This can be understood from fl3.3p by taking the hmit C{p) —>■ and taking into account the 
fact that for C{p) = 0, there is no creation so that K ^ 0. A symmetric expression holds 
for A{p) = 0. 

Note that J[o^t]{p,Q,K) = for typical evolutions, i.e. when (12.161) . (12.71) and fl2.15p are 
satisfied. 

In higher dimensions, a similar expression for the large deviation functional holds. The 
coefficients D and a should be replaced by matrices and Q becomes a vector field (see 
|BDGJLH IBDL] ). but the non-conservative part is unchanged. 

3.2. Heuristic derivation of the large deviations. We turn now to a heuristic derivation 
of the large deviation principle (13.11) . In the functional I[o,t], the contribution of the two 
currents Q{x,t),K{x,t) splits into two independent parts. 

The first term of the functional involves only the conservative part of the current and it can 
be interpreted as the local gaussian fiuctuations of Q{x, t) around its mean —D[p(x, t))dxp{x, t) 
(see (12. 7p ) with a variance a {p{x,t)) (O. Similar expressions have been obtained in 
jBDTl [BDl IBDGJLlj for diffusive type dynamics. 

The contribution $ of the non-conservative current can be understood as follows. During 
the time interval [t,t + dt], the density p{x,t) remains essentially constant in the small 
region [x, x + dx]. During [t, t + dt], particles are created in [x, x + dx] according to a process 
k~^{x) which is the sum of the independent microscopic processes at each site in the interval 
{Nx,N{x + dx)}. As the density remains essentially constant and the local equilibrium is 
maintained, the process k~^{x) is the sum of Ndx independent Poisson processes with rate 
^^c(i,?7j)(l — rji) ~ CNdx, where C = C{p{x,t)) (12.101) . Thus k'^i^x) is a Poisson process 
with a large deviation function given for large N by 

P (^j^ k+{x) = j^Ndxd?j ^ exp {-Ndxdt^Hcil^)) 



where J^* k~^{x) stands for the number of particles created in [x,x + dx] during the time 
interval [t,t + dt] and 

^c(7+)=7+logf^) - 7+ + C. (3.4) 



In the same way, particles are annihilated in [x, x + dx] according to a Poisson process 
k^{x) with rate ANdx = A[p{x,t))Ndx. When tends to infinity, the large deviations of 
k~{x) are given by the functional \1/a (13.40 . Thus observing a large deviation K{x,t) = k 
of the non- conservative current in [x, x + dx] during the time interval [t, t + dt] boils down 
to optimize the large deviations of both processes k~^{x),k~{x) under the constraint that 
k~^{x) — k~{x) = kNdxdt 

<l>{p{x, t),k)= inf {^cil'') + ^a(7")} = inf {^c(7) + ^a(7 - k)} . (3.5) 



Adding the local contributions for each space-time intervals dx and dt, we recover the non- 
conservative part of the functional (13. 2p . 
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The previous heuristics are based on the conservation of local equilibrium (thanks to the 
fast stirring) and the summation of all the local fluctuations (Gaussian for the conservative 
current and Poissonnian for the non-conservative current). A similar additivity principle has 
been already used in the framework of diffusive systems |BDlt IBD2j . A rigorous mathe- 
matical proof of the large deviations (13. ip can be found in |BL] for the microscopic model 
introduced in section [2Tn i.e. -D(p) = 1 and cr(p) = 2p(l — p) . 

3.3. Density large deviations. The large deviations for the density have been considered 
in |JLV] and these results can be recovered from the functional (13.20 . The asymptotic 
probability of observing an atypical evolution p for < t < T is given for large by 

IP[o,T] (p) ~ exp [-A^^[o,T](p)] with G[o,t]{p) = mf J[o,t](p, Q, A') (3.6) 

where the functional I[o,t] is optimized over the currents {Q, K} satisfying the conservation 
law fl216D . 

It is convenient to rewrite the conservative current in terms of a chemical potential H 

Q{x, t) = -D{p{x, t))d,p{x, t) + a {p{x, t)) d,H{x, t) (3.7) 
where H has boundary condition H{0,t) = 0. Thus (13. 2p becomes 

Ilo,T]ip,Q,K)= dt dx I ^(^(^' [d^H{x, t)f + ^ (p(x, t), k{x, t) 

Optimizing this functional over if, K under the constraint of the conservation law (I2.16p 
implies H{l,t) = at any time and 

H{x, t) = -^2$ (p(x, t), k{x, t)) , (3.8) 

where 82 stands for the derivative wrt K. From (13.50 and (13.81) . we see that there is 7 such 
that H{x,t) = log [cd^) = - log ™^ ^P^^'^^ 

k{x, t) = C{p{x, t)) exp {H{x, t)) - A{p{x, t)) exp ( - ii'(a;, t)) (3.9) 
and thus we recover the expression of |JLVj 

G[o,T] (p) = £ dt dx ^ {d^Hf + A{p) (1 - exp(-ii) {H + 1)) 

+C(p)(exp(ii)(ii-l)-l), 
where H is defined such that the following equation holds 

dtp = {D{p)d^p - a{p)d^H) + C{p) exp {H) - A{p) exp {- H) . 
The previous equation is obtained from the conservation law (I2.16P combined with (13.70 and 

dSl]). 

4. Consequences of the large deviations 

In this section, we focus on new the features of the current large deviations which do not 
exist for purely conservative dynamics. 
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4.1. Local variance. In practice large deviations are "rarely" observed, however one can 
often understand the fluctuations close to a steady state by expanding the large deviation 
functional. Using this approach, the cumulants of the current were predicted in [BD_lJ 
for general diffusive dynamics and the results were conflrmed numerically for some speciflc 
models in |BD3t IHG] . We follow the same path and focus on the second order expansion 
which is related to the variance of the current (when there is no phase transition). 

As the number of particles is non-conserved, the time integrated currents are not constant 
through the system and the deviations of the current depend on where it is measured (this 
was not the case for the conservative dynamics). In order to illustrate this fact, we are 
going to compute the deviations of the total conservative current and of a local conservative 
current (measured in a smaller region of the bulk). The mechanisms at play are different. 
In particular, to increase the current from left to right in a given region, extra particles are 
created to the left of this region and depleted to its right (see flgure[T]). This is reminiscent 
of the vortices which induce local current fluctuations in two-dimensional models |BDL] . 

Given a function A(x) with x G [0, 1], we are interested in the large deviations of the 
integrated current dx X{x)Q{x,T) = Jq dt dx X{x)Q{x,t) ovei aveTj long time interval 
[0,T] 

•^^^^ ^ T^N^ (^^P (/ dx \{x)Q{x,T)^^ . 

This boils down to minimizing the time dependent large deviation functional (13.21) 

.F(A) = hm 1 sup ( / dxX{x)Q{x,T)-I[o,T]{p,Q,K)] . (4.1) 



^ P,Q,K ^Jo ) 

For some dynamics |BD2^ IBDGJLl] . this is equivalent to a stationary variational problem 
and the supremum can be taken on time independent proflles {p,Q,K). In this case, the 
conservation law (12.161) implies K{x) = d^Q^x) and the large deviation functional (13. 2p 
reduces to 

\ iQ{x) + D(p(x))d^p(x)Y ( , , ■ , .\\ 
I{p,Q) = dx I 2a(p(x)) +$(p(x),c},.Q(x))| . (4.2) 

It is then enough to consider the simpler variational problem 

= sup I /" dx\{x)Q{x) - i{p, Q) I . (4.3) 

It is in general not true that the functional (14. ip can be reduced to (14. 3 p for any microscopic 
dynamics, i.e. for any functions Z), a. A, C and A. Some counter-examples leading to dynam- 
ical phase transitions can be found in the case of conservative dynamics in |BD2t IBDGJLl] . 
In section 14.21 we will construct a new example of phase transition in the non-conservative 
case. When the potential V (I2.12p is convex, we expect that (14. 3 p gives the correct expression 
for small deviations close to the steady state current. Indeed in this case the steady state p 
(I2.13P is unique and stable with respect to the hydrodynamic evolution (12. lip . Thus a small 
perturbation of the current should not lead to the bifurcations observed in the dynamical 
phase transitions [BD3j . 
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Remark 4.1. The previous conjecture on the equality between f l4.ip and fl4.3p can be derived 
when D is a constant, a is concave, A{p) = ap and C{p) = c(l — p) are linear (and thus 
V is quadratic). Recall that for the SSEP, one has D = 1, a{p) = 2p(l — p). We follow 
the argument in |BDGJLT] devised for the SSEP. First note that the functional I[o,t] can be 
rewritten (see e.g. |BL] ) as 



X? 



where we set for any functions f{x), q{x), k{x) with x in [0, 1] 

^(/, q) = sup ^hU, q) = sup ( / dx (g(x) + dJ{x))H{x) - \(y{f{x))H{ 
and 

$(/, k) = sup k) = sup I dx k{x)G{x) - C(/(x))(e^(^) - 1) - A(/(a;))(e-^(-') - 1) 

G G IJq 

Note that the representation holds for any A and C . 

For any H, the functional "^nif, q) is convex jointly in (/, q) because a is concave. In the 
same way, for any G, the functional k) is jointly convex in (/, k) since A and C are 

linear. As \1/ and $ are defined as a supremum of convex functions, they are also convex. 
Thus for any trajectory (p, Q, K) 

This inequality implies that the minimum in ( I4.ip can be reduced to time independent profiles. 
Indeed, any time dependent evolution (p, Q, K) which produces a total current q = -f dtQt 
has a large deviation cost greater or equal than the time independent profiles obtained by 
averaging this trajectory over time. 

Given the function A, we are going to expand T{e\) (given in (14. 3p ) to the second order 
in e. For a small shift of the current, we expect that the minimum is located close to the 
steady state p defined in ( I2.13p . We introduce the perturbation 

p(x) = p{x)+e , Q{x) = -D{p{x))d^p{x) + eq{x) , 

D[p[x)) 

with /(O) = /(I) = to satisfy the boundary conditions (12.50 . 

Expanding the functional X (14. 2 p to the second order in e leads to 

T{£\) = -ej dx \{x)D{p{x))d^p{x) (4.4) 



,./W„ \ 2a(p(x)) 2(A(p(x)) + C(p(x))) 
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where U{x) = ^ ^ by the convexity assumption of V. We introduce 

q{x) + f'{x) q'jx) + U{x) f{x) 

and optimize (I4.4p over / and q. This imphes 

X{x) = ip{x) - dr^tpix) 
= -d^^{x) + U{x)ij{x) 

with the boundary conditions ipi^) = "^(1) = and /(O) = /(I) = 0. We deduce that 

'A'(x) = U{x)iIj{x) -dli^ix)) 

(fix) = X{x) + d^ipix) (4.6) 
_V^(0) = V^(l) = /(0) = /(l) = 

Since U{x) ^ 0, the solution (v?, V') of (14. 6 p is unique. The current and density profiles 
{q, f) can then be determined uniquely since the boundary are fixed g'(0) = g'(l) = and 
/(O) = /(I) = 0. 

• Total current deviations. 

The variance of the total current Jq dxq{x) can be obtained with a constant Lagrange 
parameter A(x) = A. In this case, the solutions of (14. 6 p are (p{x) = X and ip{^) = 0. 
Plugging this in (14.40 leads to 

J'ieX) = -{eX) I dxD{p{x))d.,p{x)+^^ f dxa{p{x)) . (4.7) 

The variance of the total conservative current should be given by dxa{j){x)). 

The result is the same as the one obtained in the purely conservative case [BDl] , however 
it is important to note that the minimizers of the variational problem (14. 4p can be different 
from the conservative case and that the optimal current is, in general, no longer uniform 
through the system. 

• Local current deviations. 

In order to solve (14. 6 p for local current fluctuations, we consider the simple case where 
p is constant so that all the differential equations have constant coefficients (denoted by 
A,C,U,a). For example, this is the case if the reservoirs (12. 5p have equal densities pa = 
Pb = p and V'{p) = 0. We rewrite flTSj) and flT6l) 

'X'{x) = U-^ix) - dl (^(x)) , ^(x) = A(x) + d^il){x) , 
q{x) + fix) = a^(x), q'{x) + U f{x) = {A + C)^(x) , (4.8) 
^(0) = ^(l) = /(0) = /(l) = 0. 

Recall that by the convexity of V then U ^ 0. The solutions are given by 



tpix) = ao smh.{VUx) - dy ^ sinh [VU{x - y)) , 

(y9(x) = X{x) + aoVU cosh{VUx) — dy X'{y) cosh [Vll{x — y)) 



(4.9) 
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with „„ = r^.^^-'-Cv^C-^)). 

'q"{x) -Uq = ~Uaip{x) + {A + C)'^'{x) = {A + C - Ua)ip{x) - {A + C)X{x), 
f"{x) - U f{x) = (T^'{x) -{A + C)ij{x) = {Ua-A- C)ij{x) , 
^'(0) = g'(l) = /(0) = /(l) = 0. 
Finally the optimal current and density are given by 



q{x) = ai COsh(v^x) + /; dy iA+C-UaMy)-iA+C)X[y) gj^j^ ^^^^ _y^y 

fix) = 02 sinh(v^x) + {Ua-A-C) dy ^ sinh {VU{x - y)) , ^^'^^^ 



Vu 

with ai,a2 such that the boundary conditions are satisfied 



0.0 0.2 0.4 \ 0.6 0.8 /.O 



Figure 1. For a function A(x) (dashed line), the optimal deviation f{x) (thick 
line) from the steady p and the optimal current q{x) (thin line) have been 
rescaled to be represented on the same graph (the vertical axis unit is irrel- 
evant). The function X{x) is peaked around x = 1/2 in order to impose a 
maximum current q in the region around 1/2. The density adapts so that its 
gradient facilitates a larger flux in the region around 1/2. Particles are created 
to the left of 1/2 and removed on the right. 

When A is a constant, the optimal density and current associated to the deviations of the 
total current are constant: / = and q = Xa. However to achieve a local current deviation, 
the interplay between the stirring process and the Glauber dynamics is more complex and 
the density is no longer uniform as depicted in figure [TJ 

4.2. Phase transition. In this section, we consider a slightly different dynamics in order 
to explain in the simplest framework a mechanism of phase transition for the current large 
deviations. The system is on a periodic chain {1, N} without reservoirs. Instead of the SSEP, 
the stirring process is now the Weakly Asymmetric Exclusion Process (WASEP) with jump 
rates to the right given by A^^(l + -^)dt and to the left is given by A^^(l — -^)dt. The drift 
is of order and z/ is a fixed parameter. This stirring process is combined to the Glauber 
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dynamics with rates fl2.18l) . For 7 > 1/2, this leads to a double well potential V fl2.19p with 
minima p'^, p~ = 1 — and a local maximum at 1/2. 

The currents are defined as in section [2^ The macroscopic equation (12. 7p for the conser- 
vative current becomes 

Q{x,t) = -d^p{x,t) + ua{p{x,t)) , (4.11) 
with (j{u) = 2u{l — u). Thus the new hydrodynamic evolution can be written 

dtp{x,t) = Ap{x,t) - ud,{aipix,t))) - V'{pix,t)) . (4.12) 

The stationary measure of this dynamics exhibits a phase transition as the steady state 
concentrates on two different density profiles in the thermodynamic limit N ^ 00. To 
see this, we first notice that the constant profiles equal to p'^ , p~ and 1/2 are all stationary 
solutions of the hydrodynamic equation (I4.12p . To guess the structure of the stationary state, 
one has to consider the stochastic corrections provided by the fiuctuating hydrodynamics 
|DFLl lOSj . A noise term depending non-linearly on the local density has to be added to 
(I4.12P (see [DFLj for an exact expression). Thus only p'^,p~ are stable with respect to 
the noisy equation and 1/2 is not. Depending on the values of 7, z/ there can be also other 
stationary solutions of (I2.13P which are not constant, but they will be unstable. Thus we say 
that the stationary measure has a phase transition: the steady state concentrates in the large 

limit on the two constant profiles p"*", p~ with equal probability due to the symmetry of the 
rates fl2.18p . (The previous argument is mainly heuristic and we postpone a mathematical 
justification of the occurrence of a phase transition to a futur work.) 

We will see now that the phase transition of the stationary measure has also implications 
on the current large deviations. Adding a weak drift u to the stirring process modifies only 
the cost of the conservative current deviations in the large deviation functional (13.20 . The 
new functional reads 

, • i (Q(x,t) + d:,p(x,t) - iya(p(x,t))Y / , 
Xfo,^](p,g,K) = dtj^ dx^^:^^ 2a{p{xt)) +^[p{x^),K{x,t) 

(4.13) 

where $ is still given by (13.30 . the conductivity a{u) = 2u{l — u) and the diffusion coefficient 
D is now equal to 1. We refer to [BD2\ IBD31 IBDGJLl] for a study of the large deviations 
in conservative dynamics with weak drifts. 

As the steady states concentrates on p~^,p~, the mean conservative current (14. lip is 
z/cr(p^) = V(y{p^) (in the large A^ limit). We are going to compute the cost of an increase of 
the total current over a very large time interval. For any shift of the current q G [z/cr(p"^), z//4] 

^hm Jhn ^ logP[o,T] dxQ{x,T) ~ = . (4.14) 

The strategy to minimize the functional (I4.13P goes as follow. To produce a current close to 
i/T/4 = i>a{l/2)T over a very long time T, the system (which originally starts close to p^ 
or p^) can evolve to the constant profile 1/2 within a time of order 1. This first move has a 
cost exponential in A^ but independent of T so that it is irrelevant in the limit (I4.14p . Then 
the density remains equal to 1/2 for the rest of the time and therefore the system produces a 
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conservative current (14. Ill) asymptotically equal to z/T/4 for large T. As 1/2 is a metastable 
profile this has no large deviation cost (at least in the scaling considered in fl4.14p ). For 
an intermediate current q = aua^p'^) + (1 — a)z//4 (with a g]0, 1[), the deviation can be 
achieved by a time-dependent density profile which is close to during a time aT and then 
close to 1/2 during a time (1 — a)T. The optimal cost remains equal to 0, but the time 
dependent variational principle (I4.ip cannot be reduced to (14. 3p even for arbitrarily small 
current deviations. A flat piece in the large deviation function (see Figure [2]) is a sign of a 
phase transition. 




ua{p+) 1^/4: aT T 

Figure 2. The current large deviation function G is represented on the left. 
On the right, a switch at time aT between the densities and 1/2 is depicted. 



4.3. The fluctuation relation. The large deviation function of the entropy production 
obeys a very general symmetry property [GCl lECMj which is known to hold for a wide class 
of stochastic dynamics [Kj ILSl [M] . When applied to conservative dynamics in contact with 
two heat baths, the entropy production can be reinterpreted as the current flowing through 
the system and therefore a similar symmetry holds for the large deviation function of the 
current at the microsopic and macroscopic levels [BDlt IBDGJLl] . A generalized notion of 
detailed balance introduced in [Dl IBDSj allows to consider an arbitray number of reservoirs 
and therefore is well suited in the case of bulk creation/dissipation. Finally, the deviations 
of the entropy production for the more realistic case of granular gases was considered in 
IBGGZj . 

In this section, we recover the fluctuation relation at the macroscopic level by means of 
the functional (13.21) . 

We are going now to investigate the symmetries of the functional (13. 2p with respect to 
time reversal which is a key feature of the Gallavotti-Cohen relation. Let (p, Q, K) be an 
evolution from time to time T such that the coefficients A[p) and C{p) are positive. The 
functional (13.20 can be rewritten as 

g(x, tf+{D (p(x, t))d.,p{x, t)f ^ Qjx, t)D{p{x, t)) d^pjx, t) 
2o-(p(x,t)) cr(p(x,t)) 

+$(p(a;,t),ii:(a;,t)), (4.15) 



Xlo,T]ip,Q,K)= dt 

Jo Jo 
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with $ given by 



^p,K) = C{p)+A{p)-JK^ + AAip)Cip)~K\og 



+ AA{p)C{p) - k 



2C{p) 



where the previous identity follows from (13.31) by multiplying the numerator and the denom- 



inator (in the log) by y A'^ + AA{p)C{p) - K. 

Let S{u) be a function such that S"{u) = 2^^^ (we shall see later on that S can be 
interpreted as a relative entropy). By integration by parts, one has 

dt f dx — ^/^^/ ' l\ ' ^ = / c?^ / dxQ{x,t) dx S'(p{x,t)) 
Jo cr[p{x,t)) Jo Jo ^ 

dt dxdM^,t)S'{p{x,t)) + (^Q{l,T)S'{pt) - Q{0,T)S'{pa)) , 

where Q{y,T) = dtQ{y,t) and Pa, Pb are the densities imposed by the left and right 
reservoirs. The conservation law (I2.16P leads to 

dt I ^^^(:M)^(p(a;,t))5^'P(a;,t) 



a{p{x,t)) 

Q{l,T)S'{pb) - QiO,T)S'{pa)) + dtj^ dx{dtp{x,t)-k{x,t)^S'{p{x,t)) 

(Q{l,T)S'{pb) - Q{0,T)S'{pa)) - dt dxk{x,t)S'{p{x,t)) 
+ dx(^S{p{x,T)) -S{p{x,0)] 



Combining (I4.15p . (14.161) and the previous identity, one can compare the large deviation cost 
2^[o,T](p, Q, k) to the one of its time reversal (p, Q, k) = {p(T — t), —Q{T — t), —k(T — t)) 

j[o,T](p, Q, k) - i[o,T]{p, Q, k) = £:[o,T](p, Q^K) + dx {s{p{x, T)) - s{p{x, 0))) .(4. 17) 

where the entropy production can be identified as 

%T](p,Q,i^) = (q(1,T)5'(p,) -Q(0,T)5'(p,)) (4.18) 

^ lo L ^'(p^^''^'>) ^^^^ 



M(p(^,^)) 

\C{p{x,t)) 



By construction £^[o,t] is anti-symmetric under time reversal. 

We consider dynamics such that S is bounded (which is the case if the (microscopic) 
number of particle per site is bounded). Then the large deviation function of the entropy 
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production 



satisfies the Gallavotti-Cohen symmetry Q{e) = G{—e) — e (thanks to the identity ( 14.17P and 
the fact that the term dx (S(^p{x,T)) — S[p{x, 0)) j is bounded uniformly wrt time). 



For conservative systems, the entropy production £^[o,r] (14.181) is proportional to the con- 
servative current Q{0,T) flowing through the system and therefore the symmetry holds also 
for the large deviation functional of the conservative current |B Dll IBD31 IBDGJLl] . For 
general creation/annihilation bulk rates, the relation between ^^[o,t] and the currents is less 
straightforward as the density is coupled to the non-conservative current K. 

Simplifications occur when the steering process and the Glauber dynamics are both re- 
versible (locally) with respect to the same Gibbs measure. Let Pp be a Gibbs measure with 
constant density p which is invariant for the steering process for any density p (see (12.91) ). 
We denote by x(p) the susceptibihty of Pp. According to the fluctuation dissipation rela- 
tion [S], the function S{u) satisfies S"{p) = 2-^^ = Thus S'{p) can be interpreted as 
the chemical potential conjugated to the density p and S as an equilibrium large deviation 
function. Suppose that the creation and annihilation are reversible with respect to pp^ for a 
given density po, i.e. 

where • | {T]j}j^i) is the Gibbs measure on r]i conditionally to the configuration {r]j}j^i 
outside the site i. Thus (I2.10p can be rewritten 

Mp) = Pp{c{i, v)Vi) = exp {S'{p) - S' (po)) pp{ppo {c{i, ri)r]i \ {rij}j^i)) 
= exp{S'{p)-S'{po))Cip). 

Using this identity, the entropy production (14. ISp can be rewritten only in terms of the 
currents 

S[o,T]{p,Q,K) = Q{l,T)S'{p,) - Q{0,T)S'{pa) - S'{po) [ dxK{x,T). 

Jo 

5. Conclusion 

In this paper, we analyzed the current large deviations for non-conservative diffusive dy- 
namics driven out-off equilibrium by external reservoirs. We obtained an expression for the 
large deviation functional (13. ip which generalizes the results derived for the conservative 
dynamics |BDGJLlj - [BD3] . 

From the large deviation functional, one can understand how the deviations of the currents 
and the density are coupled. Contrary to the conservative case, there are now two currents 
instead of one and it is therefore more delicate to derive exact expressions of the minimizers 
of the functional given some constraint (for example given the total currents). It would be of 
interest to derive more systematic predictions for the cumulants of the currents as in |BDlj . 

The phase transition for the dynamics considered in this paper is currently not understood 
at the mathematical level and we plan to study it in the future. 
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In section 14.31 we showed that the entropy production fl4.18p satisfies the fiuctuation 
relation. One may wonder if other (measurable) physical quantities related to the injected 
current satisfy similar relations. In this spirit, the exact expression of the large deviation 
functional (13.11) could be helpful to solve some controversy on the validity of the fiuctuation 
relation in dissipative systems |FPj . 
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